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Abstract 

Superconducting transition temperatures (T c 's) of MgE>2, Mgi_ x Ca x .B2, 
Mgi_ x Na a .B2, and Mgi_ a: Al a: B2 are studied within the McMillan approxima- 
tion using electronic and structural information obtained from density func- 
tional theory within the generalized gradient approximation. The density of 
states and T c of MgB2 are both shown to be extremely sensitive to volume; in 
fact the density of states around the Fermi level is found to rise with increas- 
ing volume because of a prominent van Hove peak. Doping the Mg sublattice 
with small amounts of either Ca, which substantially increases the unit cell 
volume, or Na, which removes an electron from the unit cell while likewise 
increasing its volume, shifts the Fermi level toward the peak and thus both 
types of doping are predicted to enhance T c ; in Mgi_ x Al x B2, however, the 
combined effects of the additional electron and decreasing average unit cell 
volume are shown to decrease T c with increasing Al concentration, consistent 
with recent experiments. 
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I. INTRODUCTION 

Since the recent discovery of superconductivity at T c = 39 K in MgB 2 by Akimitsu!, much 
work has been devoted to uncovering the underlying superconducting mechanism. Evidence 
for conventional, phonon-mediated BCS-type superconductivity! is now abundant. A boron 
isotope effect has been observed!, and photoemission measurements indicate the opening of 
a superconducting gap in the density of states (DOS) with a narrow coherent peak!. The gap 
seems to have predominantly s-wave character!, and its temperature dependence follows the 
BCS prediction. As shown by Raman!, photoemissionB, and some tunneling experiments!, 
the ratio of the gap to T c is slightly larger than the BCS value, indicating that MgB2 is in 
the intermediate coupling regime. 

As already suggested by some authors!'! the large T c in MgB 2 can been understood 
with simple, weak-coupling arguments: the density of states at the Fermi level is large 
compared with that of an sp-valent, nearly-free electron metal of comparable density 
(r s = (SVI/Atc) 1 / 3 ~ 1.8 a.u., where Q is the volume/electron), and the phonons most strongly 
coupled to electrons around the Fermi level are due to the vibrations of the light B atoms 
alone. This simple intuition, previously used to suggest a high T c in metallic hydrogen!, indi- 
cates two routes for increasing T c : either (i) doping the B sublattice with lighter atoms or (ii) 
increasing the density of states at the Fermi level, A^(0), without significantly diminishing 
the coupling to the relevant lattice modes. Strong-coupling effects and electron-electron in- 
teractions can potentially diminish the transition temperature with respect to that expected 
from a weak coupling approximation, and therefore both require careful consideration, which 
we provide here. In this work we focus on possibility (ii). In particular, in Sec. II we first 
describe our first-principles calculations of the electronic properties of MgB2 and related 
compounds, and then discuss the details of the McMillan approximation to the Eliashberg 
equations to evaluate T c . In Sees. IIIA and IIIB we report the structural parameters, elec- 
tronic band structure, and density of states for MgB 2 as a function of volume, and for A1B 2 , 
CaB 2 , and NaB 2 at equilibrium. Importantly, a small decrease in electron density is shown 
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to increase N(0) and T c . As an example of an empirical realization of this effect, in Sec. IIIC 
we consider the result of adding small amounts of either Na, Ca, or Al to the Mg sublattice, 
which in the case of Na and Ca, decreases the electron density, or, in the case of Al, increases 
the electron density. For Mgi_ x Al x B2, we find that T c declines with increasing x in accord 
with experiments0'0. For Mg 1 _ x Na x B 2 and Mg^^Ca^E^, however, we find that T c could 
increase up to 53 K for an Na level of x = 0.2 or up to 52 K for a Ca level of x = 0.2. We 
conclude in Sec. IV. 



II. METHODS 

To investigate the structural and electronic properties of MgB2 and related com- 
pounds from first principles, we use a plane wave-implementationE! of density functional 
theoryEl within the generalized gradient approximation (GGA)0 plus Vanderbilt ultrasoft 
pseudopotentialslll. Our Mg pseudopotential treats the low-lying 2p electrons explicitly as 
valence states. For all calculations we use a dense k-point sampling (k-meshes of 19x19x15) 
and a 37 Ry plane wave cutoff, these together converging the total energy to less than 
a few meV/ion. After all lattice parameters are optimized for a given volume, electronic 
densities of states are calculated. Once the DOS is determined, we estimate T c using the 
McMillan approximation to Migdal-Eliashberg (ME) theoryi@0. The ME theory includes 
the dynamical properties of the effective electron-electron interaction (and therefore screen- 
ing of the Coulomb repulsion) by introducing the Coulomb pseudopotential /z*; self-energy 
corrections arising from the electron-phonon interaction are also included. The ME the- 
ory does not, however, include vertex corrections, since according to the Migdal theorem 



such terms are smaller than the self-energy correction by Xy'ujQ/uj e} where A is the effective 
electron-phonon coupling (A ~ 1), ujq the average phonon frequency at equilibrium volume, 
and uj e a characteristic frequency of the order of the electronic bandwidth (~ 10 eV). In the 
diboride compounds considered here, the average phonon frequency is roughly 700 K, the 
frequency of the strongly-coupled E 2g mode of MgB 2 i'0; with these numbers we estimate 
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\y/uo/u} e ~ 0.1 and therefore vertex corrections can be neglected. Moreover the DOS has 
a weak energy dependence in a range of energy of order ujq (~ 60 meV) around the Fermi 
level Ep (see Fig. 2), and hence a constant density of states N(E) = iV(0) near Ep is a 
reasonable approximation. 

With these simplifications the solution to the Eliashberg equations for T c can be ap- 
proximated by the McMillan formula (using an estimate^ for the Debye temperature, 
6 D ~ (4/3)u>„) 

T c = 0.92.„exp(-1.04 A ^ (1 + + A 062A) ), (1) 

where /i* = V c N(0)/(l+V c N(0)ln(co e /uJo)), where V c is an approximate static Fermi surface 
average of Coulomb interactions. Using u e ~ 10 eV, we obtain \n{co e /coo) ~ 5; this implies 
that V c is well renormalized, and that //* can reach a maximum value of 0.2 for large V c . 
Although a more rigorous treatment of electron-electron interactions may reduce the value of 
//* with respect to the accepted value@@, in what follows we take /i* = 0.1, the conventional 
value for sp-valent, nearly-free electron metals. Changes in V c N(0) with volume will not have 
a considerable effect on /i*, which, in these compounds, is a slowly varying function of V c N(0) 
near its limiting value. The effective electron-phonon coupling is 



_ N(0) < g 2 > 
X ~ ' (2) 



where < g 2 > is the average over the Fermi surface of the elements of the square of the 
electron-phonon matrix and M is the atomic mass, and, in this work, varies with volume 
through N(0) and uo 2 . In the following calculations of T c , the DOS is determined at each 
volume from first principles, and the volume dependence of the phonon frequency is taken to 
be u = coo(V/Vq)~ 7 , where 7 is the Gruneisen parameter and Vq is the equilibrium volume; 
the electron-phonon matrix elements and the Coulomb pseudopotential are approximated 
to be volume and concentration independent. First-principles evaluation of the volume and 
composition dependence of these matrix elements^, which will most strongly depend on 
proper treatment of the screened electron-ion interaction, is a formidable task for alloys and 
beyond the scope of this work. 



III. RESULTS 



A. Compressed and expanded MgE$2 

Under standard conditions MgB 2 crystallizes in a three atom, hexagonal primitive cell; in 
it, a single Mg atom is centered one half-lattice vector above a honeycomb boron network. 
The lattice parameters have been measured to be a = 3.084 A and c = 3.523 A at low 
temperaturesfB, and we calculate a = 3.065 A and c = 3.519 A, in excellent agreement 
with these experiments and also with previous calculations^'^. Detailed discussions of 
the electronic structure of MgB 2 have already appeared in a number of recent studies^!. 
Two dimensional hole-like cx-bands (arising from intraplanar B p x y -like orbitals) result in a 
flat DOS with a prominent van Hove spike. A three dimensional p z -\ike band also weakly 
contributes to the DOS. Our calculated bands at equilibrium volume (Fig. 1(b)) support 
these views and agree well with recent studieHHI. 

To examine the effects of volume on the electronic structure near the Fermi level, we 
optimize the structural parameters of MgB 2 at several densities, corresponding to changes 
in volume of ±30%, and calculate the energy bands and density of states; the results appear 
in Figs. 1 and 2. The cx-bands are largely nearly-free electron-like: their dispersion is 
parabolic near the V point, and their overall bandwidth is comparable to the free electron 
value at this density (~ 15.5 eV). In addition, the bandwidth increases as volume decreases, 
scaling roughly as a~ 2 as might be expected for free electrons, and the nearly-free electron 
gaps at the zone boundaries grow with increasing volume. The weak dispersion of the cx- 
bands along TA reflects their particular quasi-two dimensionality, decreasing with increasing 
volume (and increasing c-axis). 

The position and shape of these bands are reflected in the density of states appearing 
in Fig. 2, which we find to depend strongly on the cell volume. With respect to its value 
at equilibrium, N(0) decreases by as much as 24% at the highest compression considered. 
Correspondingly, as shown in Fig. 2, N(0) increases as the lattice is expanded, contrary 
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to expectations for a nearly-free electron metal. The dependence of the DOS on volume is 
almost entirely due to the changes in the width and position of a considerable van Hove 
peak (~ 2 eV below the Fermi level), which originates from a saddle point in the highest 
occupied cr-band at the M point. The decrease in bandwidth with increasing volume reduces 
the seperation between the peak and the Fermi level, driving up the density of states. The 
singularity is further enhanced by an increase in two dimensionality (due to the elongated 
c-axis evident in Fig. 2's inset), and a complementary increase in the primary Fourier 
component of the electron-ion interaction (evident from the larger gap at the zone boundary). 
Because these a-bands are essentially two dimensional, they result in a nearly constant DOS 
near the top or the bottom of the bands (where the dispersion is most free electron-like); 
their contribution to the DOS does not change appreciably with volume, and indeed, the 
DOS at each expanded volume takes on nearly the same value roughly 0.3 eV above Ep 
(Fig. 2). 

As a result of the apparent increase in N(0) for lower densities (i.e., larger volumes), 
the superconducting transition temperature T c should rise, and to verify this explicitly we 
calculate the transition temperature T c using Eqs. (1) and (2). We first determine the 
electron-phonon coupling by requiring that T c = 39.5 K at the equilibrium volume. Using 
coo = 700 K and fi* = 0.1, we find that a coupling value of A = 0.835 is requiredil, in 
agreement with recent estimatesi'illl of A between 0.7 and 0.95. We then calculate T c for 
a range of volumes (V/Vq = 0.7 — 1.3) using this value of A, with the results appearing in 
Fig. 3. A general trend is clear: at lower volumes (positive pressure) T c decreases, while as 
volume increases (negative pressure) T c increases, in accordance with the DOS at the Fermi 
level. The variation in T c does depend on the Griineisen parameter, and for larger values 
of 7, T c can change quite substantially over the volume range considered here. For example 
T c ranges from 26 to 50 K for 7 = but for 7 = 2 it can reach 79 K. A larger value of 
/i* will lower T c : our estimate of T c = 62 K at V/Vq = 1.1 with \x* =0.1 and 7 = 2 (see 
Fig. 3) drops to 51 K for [i* = 0.15 (for fixed A). As MgB 2 possesses equal numbers of 
electron and hole carriers, we also note that correlated electron-hole fluctuations may reduce 
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fi* and increase our predicted transition temperatureslll'E], as mentioned above. Hydrostatic 
pressure has already been found to decrease the transition temperature, with a slope in the 
range of -1.6 K/GPa@ and -1.1 K/GPa@. Extrapolating to low pressures from data at 
V/Vq = 0.9 (p ~ 19 GPa), we find a slope of -1.2 K/GPa for 7 = 1.5, slightly lower than a 
previous calculation!! of this slope using 7 = 1. The value of 7 is not experimentally known, 
and it is used here as an adjustable paramete ; we reserve its explicit calculation from first 
principles for a future study. 

B. Increasing or decreasing T c : the effects of doping with Ca, Na, or Al 

While the large effect of volume (and electron density) on T c is evident, isotropic expan- 
sion of the lattice is generally not possible experimentally. Epitaxial growth of MgB 2 is a 
one possible solution, and recently, Kang et a/.IH report the growth of high-quality MgB 2 
thin films on AI2O3 substrates having T c = 39 K with a sharp transition width of 0.7 K. 
Using different substrates, structures with a larger volume might be realized, though with 
different c/a ratios. A more direct method would be to replace Mg with a larger ion. Cal- 
cium, with an ionic radius^ of 0.99 A (compared to 0.65 A for Mg), would certainly increase 
the volume of the unit cell, as verified by recent first principles calculations^. On intuitive 
grounds, sodium would be a particularly attractive dopant; as suggested previously (see, e.g., 
RefsBiM!) it carries one less valence electron and would therefore add holes to the system 
(hole doping). Additionally, with an ionic radius of 0.97 A, sodium should increase the unit 
cell volume. Aluminum, on the other hand, adds electrons to the system and reduces the 
volume, and, again following intuition, small amounts of Al have already been shownffl to 
decrease T c . 

In what follows we study the behavior of T c for small amounts of Ca, Na, and Al in 
MgB2 quantitatively. Important physical changes due to the presence of these dopants are 
first elucidated through calculations of the equilibrium structural and electronic properties 
of their hypothetical bulk (MgB 2 -like) hexagonal structures. In agreement with expectations 
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based on the ionic radii and previous calculations^ 3 ^, replacing Mg with either Ca or Na 
increases the unit cell volume: our calculated volumes of Ca (36.40 A 3 ) and Na (34.11 A 3 ), 
are 27% and 19% larger than that of MgB2 (28.63) A 3 , respectively. For comparison, we 
find that AIB2 stabilizes with a primitive cell volume of 25.54 A 3 , 11% smaller than that of 
MgB 2 and in excellent agreement with the experimental value@ of 25.58 A 3 . In both NaB 2 
and AIB2 the intraplanar B-B distance decreases, but by less than 2%; in CaB2, however, 
the B-B distance increases by 5%. Adding a larger ion increases the c-axis for both Na 
(by 22% to 4.30 A) and Ca (15% to 4.069 A); in A1B 2 , however, the c-axis decreases by 
roughly 7% to 3.282 A. Thus by doping the Mg sublattice and changing the c-axis, the two 
dimensionality of the cr-bands can be either enhanced or diminished, and this will also have 
important consequences for the density of states near the Fermi level, as we now show. 

In Fig. 4(a) we plot the DOS of MgB2 and CaB2. Mg and Ca are isovalent, but according 
to the results of the previous section, adding Ca should have a large effect on the density of 
states. Indeed, CaB 2 would result in an N(0) of 0.981 states/eV-cell, an increase of 41% over 
that in bulk MgB 2 . Note the increase in magnitude of the singularity: the larger volume 
gives the cr-bands more two dimensional character, and the lower electron density results 
in a smaller bandwidth. Finally, we note that although the unoccupied cZ-states do weakly 
contribute to the density of states near the Fermi level through some hybridization with the 
p z -band, they have only a minimal effect on the cr-bands. 

In Fig. 4(b) we plot the DOS of NaB 2 , MgB 2 , and A1B 2 . Because Na, Mg, and Al are 
heterovalent, we expect changes both in the number of electrons per cell and in the volume 
to have a profound effect on the density of states near the Fermi level. In NaB2 the DOS at 
the Fermi level substantially increases, from 0.711 to 1.014 states/eV-cell; in contrast, N(0) 
drops to 0.378 states/eV-cell in AIB2. By removing one electron, the Fermi level of NaB2 
is shifted downward relative to that of MgB2 and toward the van Hove singularity arising 
from the two dimensional cr-bands. As discussed above, the larger volume again results in 
a more pronounced van Hove singularity just below the Fermi energy, as in CaB2. In AIB2 
the opposite occurs with respect to MgB 2 ; the Fermi level shifts away from a van Hove peak 
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that is now less prominent due to the increase in electron density and decrease in volume. 

The large increase in N(0) for both CaB 2 and NaB 2 clearly raises the possibility of a 
T c larger than that of MgB2. Indeed for the hypothetical CaB2 we obtain T c = 65.6 K and 
for NaB 2 we calculate T c = 68.5 K (using 7 = 0, and < g >, u , and fi* appropriate for 
MgB 2 ). But the expanded c-axis and appreciable DOS at the Fermi level, of course, call 
in to question the stability of the resulting CaB 2 or NaB 2 lattice. Neither compound, to 
our knowledge, has yet been synthesized. Therefore, instead of focusing on bulk CaB 2 or 
NaB 2 , we concentrate on the electronic properties and transition temperatures of the alloys 
M gi_ x Caa;B 2 and Mgi_a;Na x B 2 for x up to 0.2. Previous calculations^ using supercells 
constructed from the experimental lattice parameters of MgB 2 have been performed for 
large Li and Na doping fractions of x = 0.25 and 0.5. To obtain a quantitative estimate of 
T c for small x in the absence of experimental data and without resorting to large supercells, 
we calculate N(0) from first principles using a simple rigid band model that also incoporates 
the effects of the changing density due to changes in the lattice parameters. To facilitate 
comparison with Na and also test the validity of our results against existing experiments, 
we also consider Mg 1 _ a; Al a ,B 2 . 

To calculate T c , an estimate for N(0) is required for small x. As we have shown, the 
DOS of CaB 2 departs considerably from that of MgB 2 . Given that it is isovalent with Mg, 
the foremost consequence of adding small amounts of Ca for the density of states will be the 
increase in volume, though the change in electron- ion interaction will also have an effect, 
e.g., through the contribution of a small d state density near the Fermi level. In contrast, 
the most important consequence of adding small amounts of either Al or Na will certainly 
be a rigid shift of the Fermi level. But adding either Al or Na will also modify the electron 
density in two additional ways: first, as these atoms are not isovalent to Mg, the number 
of electrons will change; and second, the nominal charge of the metallic ion in the cell will 
clearly differ from that of Mg. Both will in turn affect the electron-ion interaction. The 
strong volume dependence of the DOS has already been demonstrated above; to illustrate 
the contribution of the electronic and ionic charge explicitly, in Fig. 5 we plot the band 
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structures of NaB2, MgB2, and AIB2 calculated with the lattice parameters of MgB2 (i.e., 
for fixed volume). At first glance the bands are quite similar, differing only by a shift of 
the Fermi level. Small differences are present, however, and they result from changes in the 
electron-ion interaction, the electronic density, and the ionic charge. For example, in NaB2 
the p z band crossing E F (whose constituent orbitals extend out of the B planes and closest 
to the metallic ion) shifts upward due to the weaker attraction to the ioni, requiring the 
o"-bands to shift downward relative to the Fermi level. In AIB2 the attraction is stronger and 
the p z band shifts the opposite way. Further, the gaps at the zone boundaries are largest 
for Na and smallest for Al. 

In the following determination of N(0) we explicitly include, as mentioned above, the 
effects of volume due to the dopant ion. As a working hypothesis we assume that the effects 
of the electron-ion interaction will be minimal, partially justified by the observed small 
differences in the cr-bands (which are most important to the DOS) when Mg is substituted 
for Ca at fixed volume. Furthermore, for NaB 2 and A1B 2 , we also implicitly neglect the 
effects of the extra electronic and ionic charge and perform all calculations with an Mg 
pseudopotential. In this context we note that further work on the impurity systems, using 
the virtual crystal or coherent potential approximations, would therefore be of interest. In 
order to include the effects of small dopant amounts on the average unit cell volume, we 
extrapolate linearly from the a and c determined for MgB 2 to those of CaB 2 or NaB 2 , thus 
obtaining a(x) and c(x), where x is the doping fraction. For A1B 2 , we make use of available 
experimental data, as we explain below. We then calculate the DOS of MgB 2 for a few values 
of a(x) and c(x), and, if we are considering either Al or Na, we adjust the Fermi energy by 
an amount A/i(x) so that the integral of the calculated DOS is equal to 6 + 2(1 — x) + xZj, 
where the dopant valence is 3 for Al and 1 for Na. (When treating small amounts of Ca, 
we simply take the DOS at the Fermi level of the expanded MgB 2 .) T c is then determined 
in the manner described above, keeping /i* fixed but scaling the electron-phonon coupling 
(Eq. 2) with a Griineisen parameter through the variation of uo 2 . The variation of [i* with 
volume and concentration is difficult to estimate quantitatively, particularly since although 
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-ZV(O) increases with increasing volume, the static Coulomb interaction V c declines. For these 
reasons, however, we may expect its variation with volume to be small compared with the 
electron-phonon coupling. Finally, in all calculations we approximate the electron-phonon 
matrix elements < g 2 > as fixed. 

1. M gi _ x Ca x B 2 

Table I displays the results of our calculations of the N(0) of MgB2 for lattice parameters 
corresponding to Mgi_ x Ca x B 2 up to x = 0.2. As expected, the increase in volume again 
increases iV(0), and, as can be seen from Fig. 6, this increases T c . For 7 = 2, T c can reach 
up to 52 K for x = 0.2, despite the fact that the total number of electrons and holes in 
each cell remains fixed. By calculating N(0) with Mg ions only, we neglect changes in the 
electron-ion interaction, but these should not be important to the density of states near the 
Fermi level, where the a-bands are most important. The stability of CaB 2 (i.e., x — 1) is 
of course questionable, given that Ca and B prefer to equilibrate in the cubic hexaboride 
structure^, but if it were possible to introduce even small amounts into the MgB 2 lattice, 
an appreciable rise in T c could be detected. 

2. M gi _ x Na x B 2 

In Table II we report DOS and lattice parameters as a function of Na doping fraction 
x at the adjusted Fermi level of MgB 2 , A/i(x). Note that N(Afi(x)) increases by 7% for 
x = 0.2 with respect to x = 0. This conflicts with the calculations of Medvedeva et alM, 
who report no increase in DOS at the Fermi level at x = 0.25 and only a 4% increase at 
x = 0.5 in their supercell calculations without relaxation of the lattice parameters. Failure 
to treat the effects on the DOS due to the increase in volume is the likely reason for this 
discrepancy. In Fig. 7 we plot T c as a function of Na doping fraction x for different 
Griineisen parameters. For 7 = 2, a T c of ~53 K would be expected for 20% Na. Although 
we explicitly considered the important effects of changing volume on N(0), as mentioned 
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above we neglect the effects of the change in ionic charge and the contribution of the reduced 
number of electrons to the electron density. Conveniently these effects should be small and 
may even largely counterbalance one another, as follows. Replacing the Mg ion with Na 
will shift the average position of the a- and p z -bands with respect to the Fermi level, as 
discussed above, but, as we have emphasized, N(0) is most sensitive to the van Hove peak 
arising from the cx-bands. Although this electrostatic effect shifts the Fermi level slightly 
away from the peak in the case of Na doping, the larger van Hove peak due to the increased 
two dimensional character should compensate. Further experiments will play a crucial role 
in determining the feasibility of Na doping, but our work nonetheless serves to illustrate 
the positive effects of increasing volume on the transition temperature, even for small Na 
concentrations. 



3. M gi _ x Al x B 2 

To examine the effects of small amounts of Al on T c we use the lattice parameters and 
electron density measured by Bianconi et a/.0; we then calculate the density of states at 
the Fermi level and T c as above!!, using an Mg pseudopotential. From Table III it is clear 
that doping the Mg sublattice with Al decreases the DOS at the Fermi level, and as is 
apparent from Fig. 8, the T c 's we calculate reproduce the data quite well using a Griineisen 
parameter of 1.2. Larger (smaller) values of 7 shift the curves to lower (higher) temperatures 
but do not change the overall shape, which is dominated by the changes in the density of 
states. Simply shifting the Fermi level without changing the volume (and therefore with 
7 = 0) results in a T c that overestimates the data considerably at larger concentrations (by 
more than 3 K at x = 0.2). In the case of Al doping the secondary electrostatic effect due 
to the increase in ionic charge would shift the Fermi level toward the van Hove peak, but 
as the doping fraction increases, the position of this peak is less relevant, since electrons 
first fill the <7-bands, shifting the Fermi level into an energy range dominated by the three 
dimensional p z band, whose contribution to the DOS is significantly less than that of the 
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2D cr-bands. Thus the contribution from the changing DOS (via rigid shift and changing 
volume) is evidently most important to T c for Al concentrations up to 8%. The additional 
electrostatic consequences stemming from the increase in ionic charge and electron density 
are artificially absorbed into 7, but we expect these effects to be minimal for the reasons 
argued above. 



IV. CONCLUSIONS 

We have found here that the superconducting transition temperature is quite sensitive 
to the unit cell volume and in particular, increasing the volume should increase T c . For the 
isotropically expanded lattice, the van Hove peak from the two dimensional bands increases 
the density of states around the Fermi level. The importance of this effect has also been 
demonstrated through calculations of three other sp-valent diborides, CaB 2 , NaB 2 , and A1B 2 , 
and also through the alloys Mg 1 _ 2 .Ca a; B2, Mg 1 _ a .Na :r B2, and Mg1_a.ALt.B2. Additional small 
amounts of Ca are found to increase T c up to 52 K due to the increase in volume without 
valence change. Strontium or barium, other possible dopants likely to increase the volume, 
would also be of interest, though low-lying d states may play a more significant role near 
the Fermi level. The presence of Na impurities in small amounts is also shown to increase 
N(0) and T c up to 53 K, contrary to the theory of hole superconductivity^. Lithium is 
another possible hole dopant, but since we calculate the equilibrium volume of LiB 2 to be 
8.9% less than that of MgB 2 , the decrease in volume may compensate for the shift in the 
Fermi level toward higher density of states, resulting in only a small variation in T c . Doping 
with Al is found to decrease N(0) and T c for small x, and a rigid shift of the Fermi level of 
a compressed MgB 2 lattice results in good agreement with recent experiments. Finally, we 
hope that this work will motivate experimental studies of impurities in MgB 2 that decrease 
the average electron density. 
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TABLES 



Table I. Lattice parameters a and c/a, and density of states at the Fermi level iV(0) 
(states/eV-cell) for different Ca concentrations x in Mgi-^Ca^L^ as calculated in this work. 



x a (A) c/a N{0) 

0.00 3.065 1.148 0.711 

0.05 3.072 1.154 0.714 

0.10 3.080 1.160 0.717 

0.15 3.087 1.167 0.723 

0.20 3.095 1.173 0.727 



Table II. Lattice parameters a and c/a, total number of electrons/cell N e , Fermi level shift 
A/x, and density of states at the shifted Fermi level N(A/x) (states/eV-cell) for the different 
Na concentrations x in Mgi^Na^E^ as calculated in this work. 

x a (A) c/a N e N(A/x) (eV) 

0.00 3.065 1.148 8.00 0.711 0.000 

0.05 3.063 1.162 7.95 0.724 0.069 

0.10 3.061 1.175 7.90 0.739 0.139 

0.20 3.057 1.203 7.80 0.760 0.282 
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Table III. Lattice parameters a and c/a, total number of electrons/cell N e , Fermi level shift 
A/z, and density of states at the shifted Fermi level N(A/x) (states/eV-cell) for different Al 
concentrations x in Mgi_ x AL E B2. The x, a, c/a, and N e are taken from Bianconi et a/.0; in 
particular note that we use the experimental lattice parameters even in the case of x = to 
better compare with experiment. N(A/i) and A/i are calculated as described in the text. 



x a (A) c/a N e N(A/x) (eV) 

0.00 3.085 1.141 8.00 0.720 0.000 

0.02 3.084 1.140 8.02 0.716 0.028 

0.04 3.083 1.139 8.04 0.709 0.063 

0.08 3.081 1.135 8.08 0.700 0.115 
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FIG. 1. Energy bands of MgB2 (a) after expansion (V/Vo = 1.3), (b) at equilibrium (V/Vo = 
1.0), and (c) under compression (V/Vo = 0.7). Vo is the equilibrium volume, which we calculate 
to be 28.63 A 3 . 
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FIG. 2. Density of states of MgB2 as a function of volume; Vq = 28.63 A 3 is the calculated 
equilibrium volume.) The inset shows the calculated c/a ratio as a function of volume. At each 
volume both a and c/a are fully relaxed. We obtain a = 2.797 A and c/a = 1.0579 for V/Vq = 0.7 
~ 97 GPa); and we calculate a = 3.284 A and c/a = 1.213 for V/V Q = 1.3. 
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FIG. 3. Calculated transition temperatures determined from the McMillan equation using 
fj,* = 0.1 as a function of volume for MgB2 with different Griineisen parameters 7. 
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E - E F (eV) 

FIG. 4. Density of states of (a) MgB2 and CaB2 at their calculated equilibrium volumes. 
Density of states of (b) NaB2, MgB2, and AIB2, also at their calculated equilibrium volumes. 
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FIG. 5. Bands of (a) NaB2, (b) MgB2, and (c) AIB2 calculated with the structural parameters 
of MgB2- With the volume fixed, the most noticable difference between bands is a rigid shift 
of the Fermi level, though additional dissimilarities are evident from changes in the electron-ion 
interaction, electron density, and ionic charge as Mg is replaced with Na or Al. 
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FIG. 6. Predicted transition temperatures (calculated with the McMillan equation using 
= 0.1) as a function of Ca concentration x for Mgi_ x Ca x B2 with different Griineisen parameters 
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FIG. 7. Predicted transition temperatures (calculated with the McMillan equation using 
= 0.1) as a function of Na concentration x for Mgi_ x Na x B2 with different Griineisen parameters 
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FIG. 8. Comparison of measured transition temperatures of Mgi_ a; Al a ;B2 as a function of Al 
concentration (x) from RefJHI with those calculated from DOS obtained from a rigid band model 
(RBM) (open squares) and from the parameters in Table III with corrections for the increase in 
volume ("RBM+volume correction") (open triangles). The lines simply serve to guide the eye. 
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